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nonlinear  response.  The  measurements  also  included  input  current,  resonance 
frequency,  and  mechanical  quality  factor  versus  electric  driving  field. 


A nonlinear  theory  for  interpretation  of  the  measured  data  was  developed  for 
the  piezoelectric  ring,  and  the  resulting  differential  equation  of  motion  was  programmed 
in  an  analog  computer.  The  output  plots  showed  nonlinear  •jump*  phenomena  in  the 
frequency  region  about  resonance  as  well  as  *beat*  phenomena  in  the  current  envelope 
response,  both  of  which  were  observed  experimentally.  Thereto  re,  analog  computer 
techniques  may  be  the  most  convenient  means  of  interpreting  and  understanding  the 
nonlinear  behavior  of  ceramic  transducer  materials. 
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DYNAMIC  NONLINEARITIES  IN  PIEZOELECTRIC  SONAR  CERAMICS 

INTRODUCTION 

Transducers  for  Navy  sonar  systems  are  often  operated  at  such  high 
power  levels  that  the  polarized  ferroelectric  ceramics  in  them  begin  to 
exhibit  nonlinear  behavior.  The  degree  of  nonlinearity  depends  on  such 
factors  as  electric  field  level,  dynamic  mechanical  stress  level,  and 
static  bias  stress  as  well  as  on  temperature  and  time.  The  overall 
influence  on  the  ceramic  in  a practical  transducer  is  a combination  of 
these  factors  and,  as  such,  presents  a complex  situation  for  interpre- 
tation of  nonlinear  behavior.  To  alleviate  this  complex  situation, 
experiments  are  usually  designed  to  treat  the  various  influential  fac- 
tors independently,  and  measurement  procedures  and  representative  ex- 
perimental data  for  ferroelectric  nonlinearities  in  ceramics  can  be 
found  in  several  references.1’2 

If  transducer  ceramics  are  to  be  operated  at  their  ultimate  power 
capacities,  substantial  nonlinear  effects  will  have  to  be  tolerated  and 
a nonlinear  design  theory  will  be  needed.  The  alternative  is  to  reduce 
the  nonlinear  effects  in  currently  used  ceramics  (for  example,  by  pre- 
stressing) or  to  develop  new  improved  ceramic  materials  that  are  linear 
over  a wider  range  of  transducer  operating  conditions.  In  either  case, 
it  is  essential  to  develop  a better  understanding  --  through  experimen- 
tation and  theoretical  analysis  --  of  the  nonlinear  behavior  of  ceramics 
and  to  somehow  relate  the  findings  to  the  microscopic  structure  of  the 
ceramic  materials. 

The  investigation  of  nonlinear  behavior  in  piezoelectric  ceramics 
has  been,  for  many  years,  an  intriguing  area  of  research  but  without 
fruition.  Even  with  the  use  of  acoustic  dummy-load  techniques,  where 
practical  design  data  for  transducer  ceramics  may  be  obtained  under 
simultaneous  application  of  high  electric  field  and  high  stress,  both 
static  and  dynamic,  the  basic  underlying  mechanisms  responsible  for  the 
nonlinear  behavior  are  not  well  understood.  For  example,  a nonlinear 
"jump"  phenomenon  has  been  detected3  in  ferroelectric  ceramics,  but 
methods  for  explaining  the  behavior  are  nil  even  though  the  behavior  of 
the  nonlinear  oscillator  is  well  known.4  This  disparity  is  attributable 
to  the  fact  that  very  little  simplified  nonlinear  theory  applicable  to 
the  low  frequency  domain  (as  opposed  to  high  frequency  nonlinear  theory 
that  can  be  validated  by  ultrasonic  measurements)  and  of  practical  util- 
ity for  transducer  ceramics  has  been  adapted  to  the  ferroelectric  ce- 
ramic materials.  A simplified  mathematical  approach  is  preferred  as  a 
tool  at  present  because  of  the  complexity  that  would  exist  in  analyzing 
large  transducer  arrays  if  a continuum-physics  approach5  were  used. 
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Therefore,  the  initial  efforts  in  this  study  were  directed  at 
developing  a dynamic  nonlinear  theory  for  a simple  structure  --  a thin- 
walled,  axially  (longitudinally)  polarized  ceramic  ring  that  undergoes 
essentially  a single-degree-of-fxeedom  vibration.  This  particular  shape 
was  chosen  for  convenience  of  measurement  as  well  as  for  ease  of  apply- 
ing nonlinear  theory. 

'•.esults  of  the  investigation  are  presented  in  terms  of  measurement 
procedure,  representative  data,  experimental  data  and  interpretation, 
and  nonlinear  theoretical  analysis. 

MEASUREMENT  PROCEDURE 

The  intent  was  to  provide  a measurement  system  capable  of  deter- 
mining simultaneously  as  many  dynamic  variables  as  possible  for  piezo- 
electric ceramic  specimens  driven  electrically  at  their  fundamental 
resonance  frequencies  in  a "free"  condition  (minimum  mechanical  damping) 
in  discrete  steps  from  low-level  values  up  to  and  including  the  value 
at  which  fracture  of  the  specimens  occurs. 

SAMPLES 

The  piezoelectric  material  initially  chosen  for  investigation  was 
Navy  Type- 1 6 ceramic,  a hard  lead-zirconate-titanate  composition,  be- 
cause of  its  extensive  usage  in  present-day,  high  power  sonar  trans- 
ducers. The  ceramic  specimens  were  in  the  form  of  thin-walled  rings 
with  dimensions 


Outside  diameter 

2.50 

in. 

(63.5 

urn) 

Wall  thickness 

0.25 

in. 

(6.35 

Dm) 

Height 

0.50 

in. 

(12.7 

Dm) 

Ratio  of  height  to 

mean  radius 

0.444 

Ratio  of  outside  diameter 

to  wall  thickness 

10. 

The  dimensions  are  sufficient7  to  ensure  pure  hoop-mode  resonance 
of  the  ring,  and,  therefore,  the  rings  may  be  considered  to  be  under- 
going a single-degree-of-freedom  vibration. 

The  ring  samples  measured  in  this  investigation  had  previously 
been  measured8  to  determine  compliance  with  the  standard  low-level 
parameter  values  for  the  material.6 

ELECTRICAL  CONNECTIONS 

The  rings  have  fired-on  silver  electrodes  on  the  upper  and  lower 
surfaces  and  are  axially  polarized.  The  driving  voltage  is  applied  to 
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the  rings  through  thin  (0.005  in.  (127  pm)  diameter)  single-strand 
wires.  The  connections  were  initially  made  by  simply  soldering  the 
wire  leads  to  the  silver  electrodes,  but  that  arrangement  appeared  to 
cause  localized  heating  effects  (two  trial  rings  fractured  approximate- 
ly at  the  solder  joints)  and  to  lead  to  separation  of  the  connections 
at  large  displacement  amplitudes.  The  final  method  of  making  the  con- 
nections therefore  consisted  of  soldering  the  electrical  lead  to  a small 
copper  tab,  cementing  the  tab  to  either  the  upper  flat  electrode  (high 
electric  potential  side)  or  the  curved  outside  surface  of  the  ring  near 
the  bottom  electrode  (low  electric  potential  side),  and  ensuring  elec- 
trical continuity  by  applying  air-dry  silver  paint  between  the  tab  and 
the  silver  electrode  surface.  (A  tab  or  electrical  lead  could  not  be 
connected  directly  to  the  bottom  electrode  surface  of  the  ring  because 
there  was  insufficient  room  between  the  surface  and  the  aluminum  base 
of  the  gas-bearing  suspension  system,  which  also  served  as  electrical 
ground . ) 

SUSPENSION 

The  rings  were  mounted  freely  on  a gas-bearing  suspension  system^* 
(see  figure  1)  to  minimize  mechanical  damping  during  the  dynamic  meas- 
urements; the  assembly  is  specifically  designed  for  piezoelectric  ce- 
ramic rings  with  an  inside  diameter  slightly  larger  than  the  central 
support.  Vertical  levitation  is  provided  by  a precision  circular  gas 
channel  in  the  aluminum  bearing  base,  and  horizontal  stability  is  pro- 
vided by  a circumferential  gas  channel  located  in  the  nylon  cylindrical 
support  on  top  of  the  base. 

Each  channel  was  fed  from  its  own  gas  reservoir  and  pressure  was 
applied  to  the  rings  via  three  small  orifices,  about  0.008  to  0.010  in. 
(200  to  250  ym)  in  diameter  and  spaced  12C°  apart  in  each  channel.  Two 
separate  regulators  (fed  from  a single  nitrogen  supply  tank)  controlled 
the  pressure  in  the  gas  reservoirs  and  were  set  at  7 lbf/in.2  and  9 lbf/ 
in.2  (48  and  62  kN/m2)  for  the  levitation  and  stabilizing  channels, 
respectively;  these  settings  produced  the  best  ring  stability. 

The  vertical  suspension  of  a single  ring  (weight  approximately  0.1 
kg)  was  approximately  0.0005  to  0.0010  in.  (12.7  to  25.4  urn)  above  the 
aluminum  base.  This  was  the  reason  the  lower  electrical  lead  was  placed 
on  the  side  of  the  ring  and  not  directly  on  the  lower  silver  electrode. 
The  gap  spacing  between  the  inside  surface  of  the  ring  and  the  nylon 
support  ranged  from  0.005  to  0.010  in.  (127  to  254  pm),  depending  on  the 
degree  of  eccentricity  of  the  ceramic  ring. 
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Figure  1.  Gas -Bearing  Suspension  System 


DISPLACEMENT  MEASUREMENT 

The  measurement  system  is  based  on  an  instrument  called  the  Fotonic 
Sensor,  manufactured  by  Mechanical  Technology,  Inc.  The  light-probe 
assembly  can  be  seen  behind  the  ring  in  figure  1,  and  the  displacement 
measurement  apparatus  is  shown  completely  in  figure  2.  The  probe,  dia- 
gramed in  figure  3,  is  composed  of  two  bundles  of  optical  fibers  --  one 
to  transmit  light  to  a reflecting  surface,  and  the  other  to  return  re- 
flected light  to  a photodetector.  The  fibers  at  the  end  of  the  probe 
are  randomly  distributed  inside  a metallic  cylinder  approximately  3.0 
in.  (76.2  ran)  long  by  0.125  in.  (3.18  ran)  in  diameter.  The  face  of  the 
cylinder  containing  the  random  distribution  of  fibers  is  positioned  (by 
means  of  a rack  and  pinion  mount)  near  the  surface  to  be  measured,  and 
the  separation  between  the  face  of  the  probe  and  the  surface  is  deter- 
mined with  a calibrated  dial  indicator  (see  figure  5).  Thus,  measure- 
ments can  be  made  without  contacting  or  mechanically  loading  the  object. 
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ZONE  OF 
ILLUMINATION 


Figure  3.  Diagram  of  Light-Probe  Assembly 


Figure  4.  Typical  Calibration  Curve  for  Fotonic  Sensor 
(sensitivity  (mV/nm)  = tan  6 ; 1 mil  = 25.4*im) 
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Calibration  of  the  instrument  is  done  statically,  and  a typical 
curve  is  shown  in  figure  4.  A large  linear  region  of  operation  is 
available,  and  for  vibration  measurements  the  distance  of  the  optical 
probe  from  the  surface  of  the  object  is  approximately  the  center  of  this 
linear  region.  The  slope  of  the  calibration  curve  is  the  dynamic  sensi- 
tivity of  the  instrument  and  is  dependent  on  the  reflectivity  of  the 
surface  of  the  object.  Because  the  ceramic  itself  is  not  a good  reflec- 
tor, a small  polished  copper  tab  was  placed  on  the  outside  cylindrical 
surface  of  each  ring  to  enhance  the  dynamic  sensitivity. 

The  particular  model  of  sensor  used  in  these  measurements  had  a 
flat  frequency  response  from  dc  to  100  kHz.  Figure  5 is  a diagram  of 
the  ceramic  ring,  the  gas-bearing  assembly,  and  the  displacement  meas- 
urement equipment.  / pf 
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Figure  5.  Block  Diagram  of  Dynamic-Displacement 
Measurement  Apparatus 


TR  5263 


ELECTRICAL  DRIVING  SYSTEM 

Figure  2 shows  the  complete  measurement  system.  The  large  unit  in 
the  top  of  the  rack  is  a complex  impedance-admittance  meter  (CIAM)  man- 
ufactured by  Dranetz  Engineering  Laboratories,  Inc.;  it  is  used  to  meas- 
ure the  low-level  linear  parameters  of  the  ceramic  test  specimens.  This 
unit  operates  with  a constant  continuous  wave  (CW)  signal  amplitude  of 
3.0  Vrms  across  the  ceramic.  Low-level  pulse  measurements  at  the  same 
voltage  are  compared  with  the  CW  measurements  to  ensure  that  the  pulse 
system  is  operating  properly.  Figure  6 is  a block  diagram  of  the  elec- 
trical driving  system,  which  utilizes  three  separate  power  amplifiers 
in  order  to  provide  clean  electrical  driving  signals  from  low  level  to 
high  level  (fracture) . A sinusoidal  driving  voltage  of  short  duration 
(pulsed  signal)  is  used  to  keep  heating  effects  in  the  ceramic  to  a min- 
imum, especially  at  high  levels.  A Tektronix  549  storage  oscilloscope 
retains  the  voltage,  current,  and  displacement  pulse  envelopes  and  the 
expanded  waveforms,  and  a Polaroid  oscilloscope  camera  is  used  to  perma- 
nently record  the  information  for  analysis.  A coherent  decade  frequency 
synthesizer,  type  1161-A  manufactured  by  the  General  Radio  Company,  is 
the  signal  source  because  of  its  frequency  stability  and  low  harmonic 
content  when  high  mechanical  quality  factor  materials  are  being  investi- 
gated . 
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Figure  6.  Block  Diagram  of  Electrical  Driving  System 


TR  5263 


REPRESENTATIVE  DATA 


Some  of  the  basic  types  of  information  available  with  the  present 
measurement  setup  are  illustrated  in  figure  7.  Figure  7a  shows  (from 
top  to  bottom)  the  profiles  of  voltage,  current,  and  displacement  pulses 
for  a trial  run  at  low  drive  level  --  approximately  3.5  V,  or  7 x 10" 3 
V/mil  (0.28  kV/m)  across  a 0.50  in.  (12.7  ram)  ceramic  specimen.  The 
time  scale  on  the  horizontal  sweep  of  the  oscilloscope  is  20  msec  per 
division,  and  the  voltage  pulse  length  shown  is  approximately  80  msec. 
Close  scrutiny  of  the  lower  displacement  trace  indicates  that  steady- 
state  conditions  exist  over  the  range  from  50  to  80  msec.  In  other 
words,  a minimum  driving  pulse  length  of  50  msec  is  required  for  the 
ring  to  reach  steady-state  conditions,  and  driving  it  with  a voltage 
pulse  less  than  50  msec  long  will  not  yield  information  indicative  of 
steady-state  material  parameters,  --  which  is  one  goal  of  this  investi- 
gation. 


Figure  7b  shows  an  expanded  view  of  the  current  decay  after  voltage 
cutoff;  the  current  was  determined  from  the  voltage  drop  across  a 1-ft 
resistor  in  series  with  the  ring  impedance  at  resonance  (see  figure  6; 
the  impedance  of  the  ring  is  roughly  100  times  greater  than  that  of  the 
1-fi  series  resistor),  from  which  the  mechanical  quality  factor  Qm  of  the 
ring  was  obtained.  The  current  decay,  which  is  indicative  of  the  me- 
chanical losses  in  the  ceramic,  is  not  strictly  exponential  at  high 
drive  levels  (as  is  common  in  linear  systems)  and  offers  the  greatest 
potential  for  measurement  error  and  misinterpretation  in  the  overall 
setup.  Gerson10  employed  a resistance  substitution  method  to  demon- 
strate to  some  extent  that  Qjq  decreased  with  increasing  stress  level 
and  that  the  stress  level  depended  on  the  electrical  drive  level. 
Therefore,  when  the  voltage  is  abruptly  shut  off  at  the  end  of  the  pulse, 
the  stresses  in  the  ring  continuously  decrease,  causing  the  so-called 
quality  factor*  to  apparently  increase  with  time  of  decay.  Although  a 
more  detailed  discussion  of  this  specific  problem  will  follow,  let  it 
suffice  at  this  point  to  say  that  the  decay  of  the  current  after  voltage 
shutoff  in  a steady-state  condition  is  normally  used  as  a measure  of 
the  losses  in  a linearly  responsive  system. 


♦Quality  factor  is  defined  in  electrical  engineering  terminology 
for  linear  systems  wily  and,  as  such,  is  undefined  when  nonlinear 
phenomena  exist. 
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Figure  7c  shows  the  steady-state  magnitudes  of  all  three  quantities 
(voltage,  current,  and  displacement)  at  resonance  of  the  ring  by  time- 
delaying  a 50  ysec/cm  portion  of  the  steady-state  region  of  the  pulse 
envelope.  These  magnitudes  are  used  to  plot  current  and  strain  values, 
as  well  as  stress  levels,  versus  electric  driving  fields,  irrespective 
of  phase  relationships. 

Figure  7d  depicts  the  magnitudes  of  the  voltage  and  current  at  1 
kliz,  a frequency  that  is  well  below  resonance  of  the  ring  (which  occurs 
at  approximately  18.64  kHz)  in  order  to  ensure  being  in  the  stiffness- 
controlled  frequency  region  to  obtain  capacitance  values  at  discrete 
drive  levels.  This  measurement  yields  the  variation  of  the  permittivity 
with  drive  level  and  can  be  used  to  correct  the  calculated  stress  level 
in  the  ring. 

Additional  information  of  the  type  shown  in  figures  7e,  7f,  and  7g 
is  also  available.  The  figures  are  traces,  with  the  same  vertical  scal- 
ing previously  used,  of  the  current  through  an  Edo  Western  ceramic  ring 
at  a driving  field  of  0.1  V/mil  (3.93  kV/m) . Figure  7e  is  the  current 
at  the  resonance  frequency  f0  and  shows  that  the  initial  current  res- 
ponse to  this  field  level  does  not  have  the  same  uniform  characteristics 
as  the  response  at  low  level  (see  figure  7a);  the  resonance  frequency 
has  also  decreased  about  150  Hz  from  the  initial  low-level  value  (18.6 
kHz)  for  this  ring.  Figures  7f  and  7g  show  the  current  responses  at 
frequencies  fj  and  fu,  10  Hz  below  and  above  resonance,  respectively. 
These  traces  demonstrate  the  highly  nonlinear  behavior  of  Type- I ceramic 
at  large  vibration  amplitudes  inasmuch  as  the  response  is  not  symmetri- 
cal in  the  frequency  domain  about  resonance.  As  the  amplitude  is 
increased  further,  jump  phenomena,  characteristic  of  a nonlinear  oscil- 
lator,11 occur,  as  demonstrated  by  Negishi.3  The  mechanical  displace- 
ment amplitude  of  the  ring  exhibits  the  same  phenomenon. 

Figure  8 illustrates  the  benefits  of  filtering.  The  top  trace  in 
figure  8a  is  the  output  of  the  Fotonic  sensor  and  shows  the  low  frequen- 
cy distortion  caused  by  erratic  motion  of  the  electrical  leads  attached 
to  the  ring  and  of  the  ring  itself  on  the  gas  bearing.  The  bottom  trace 
is  the  same  output  after  filtering  through  a 2600  Hz  to  50  kHz  bandpass 
filter;  figure  8b  is  an  expanded  view  of  the  steady-state  signal  and 
shows  the  amount  of  attenuation  introduced  by  the  filter. 


EXPERIMENTAL  DATA  AND  INTERPRETATIONS 

Several  rings  were  driven  electrically  to  dynamic  fracture,  but 
only  two  rings  (Edo  Western  ring  number  001  and  Channel  Industries  ring 
number  23,  henceforth  referred  to  as  ring  A and  ring  B,  respectively) 
were  measured  sufficiently  to  provide  reliable  data. 


11 


TR  5263 


8a.  Unfiltered  and 
filtered  envelopes 


8b.  Unfiltered  and 
filtered  steady-states 


Figure  8.  Fotonic  Sensor  Output  Waveforms 


The  measurements  were  recorded  under  pulsed  electrical  drive  with 
the  rings  mounted  on  the  gas-bearing  suspension  system.  Initially,  low- 
level  pulse  measurements  were  made,  and  the  results  were  compared  with 
continuous  wave  measurements  made  at  the  same  drive  level  in  order  to 
ensure  that  the  pulse  measurement  system  was  operating  properly.  Dif- 
ferences in  the  resonance  frequencies  for  the  two  methods  were  less  than 
0.1  percent,  and  the  values  of  Qm  determined  from  the  electrical  admit- 
tance plot  and  from  the  current  decay  were  within  2 percent  of  each 
other.  This  latter  result  gave  credence  to  the  technique  of  determining 
Qm  from  the  current  decay.  The  measured  mechanical  displacement  of 
ring  A was  within  10  percent  of  the  value  predicted  from  the  low-level 
linear  parameters  for  that  specimen  and  verified  the  calibration  of  the 
displacement  probe. 

SUSPENSION-SYSTEM  LOSS 


During  this  low-drive-level  phase  of  the  investigation  it  was  dis- 
covered that  the  gas-bearing  suspension  system  introduced  unaccountable 
mechanical  losses.  This  problem  is  evidenced  by  figure  9,  which  shows 
two  complex  admittance  plots  of  ring  B — one  with  the  ring  freely  sus- 
pended from  a rubber  band,  and  the  other  with  the  ring  mounted  on  the 
gas  bearing.  Under  the  low-level  linear  drive  conditions  used  for  the 
measurements,  Qm  is  defined  as 


Qm  = 


u)E  M 
r 


f2-fl 


(1) 
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where 


is  the  resonance  frequency  at  constant  electric  field  E, 

M is  the  mass  of  the  ceramic  ring  (M  = p V,  where  p is  the 
density  of  the  material  and  V is  the  volume  of  the  ring), 

R^  is  the  mechanical  loss  associated  with  the  ring  material,  and 

fj,  f^  are  the  quadrantal  or  half-power  frequencies. 


f6 

r 


Figure  9.  Complex  Adnittance  Plots  for  Hard  Lead-Zirconate-Titanate 

C 

Ceramic  Ring  (f£  is  hoop-mode  resonance  frequency; 
fl  and  f 2 are  quadrantal  frequencies) 
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When  the  ring  was  mounted  on  the  gas  bearing,  the  value  for  de- 
creased 23  percent  (that  is,  the  electromechanical  resistance  increased 
from  97.5  R to  130.7  R) , and  the  resonance  frequency  increased  about 
15  Hz.  These  anomalies  are  probably  due  to  the  fact  that  there  is  a 
confined  air  space  between  the  inside  wall  of  the  ring  and  the  nylon 
center  support  of  the  gas  bearing  that  introduces  acoustic  stiffness 
and  loss.  Improvements  in  the  design  of  the  gas  bearing  should  allevi- 
ate such  problems.  Inasmuch  as  there  was  no  apparent  correction  factor 
that  would  account  for  the  additional  loss  and  inasmuch  as  the  gas- 
bearing suspension  was  still  needed  to  make  displacement  measurements 
on  the  rings,  all  variations  in  parameters  at  higher  levels  were  com- 
pared with  low-level  values  measured  with  the  rings  mounted  on  the  same 
gas  bearing. 


OBSERVED  PARAMETERS 


Table  1 lists  the  raw  data  measured  for  rings  A and  B.  The  volt- 
age, current,  mechanical  displacement,  and  resonance  frequency  were 
measured  directly,  and  Qm  was  determined  from  the  decay  of  the  current 
envelope.  An  accurate  determination'  of  Qm  is  hampered  by  the  fact  that 
once  the  driving  voltage  is  terminated,  the  amplitude  of  vibration 
(which  is  proportional  to  the  stress  and  strain  in  the  ceramic  ring) 
does  not  decrease  in  a simple  exponential  manner  and  therefore  causes 
the  apparent  Qm  to  vary  with  time  of  decay.  Figure  10  demonstrates 
emphatically  thq.  degree  of  variation  experienced  on  ring  A at  an  elec- 
tric field  level  of  0.22  V/mil  (8.7  kV/m) . The  apparent  Qm  increased 
approximately  180  percent  within  a time  span  of  20  msec.  Therefore,  in 
order  to  alleviate  the  situation  and  to  obtain  a value  considered  to  be 
representative  of  the  actual  losses  of  the  ceramic  under  the  steady- 
state  driving  conditions  at  the  designated  levels,  Qm  values  for  all 
drive  levels  were  determined  from  the  first  few  points  on  the  decay 


For  future  measurements,  the  accuracy  can  be  improved  by  expanding 
the  region  of  the  decay  curve  near  voltage  cutoff  and  measuring  the 
decay  over  the  first  few  cycles.  In  addition,  a resistance  substitu- 
tion method10  could  be  employed  during  the  steady-state  portion  of  the 
driving  pulse  in  order  to  determine  the  relative  differences  between 
the  two  methods.  These  approaches  were  not  utilized  for  this  series  of 
measurements  because  the  variation  problem  in  Qm  with  decay  at  high 
drive  level  was  not  recognized  until  after  the  measurements  had  been 
completed  (and  the  samples  had  fractured). 
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Table  1.  Measured  Data  for  Ceramic  Rings 


Ampli- 

fier* 

Electric 

Field 

(kV/m) 

Current 
I (mA) 

Mechanical 

Displacement 

ur=b(um) 

Resonance 

Frequency 

f^  (kHz) 

Mechanical 

Quality 

Factor 

Qm 

0.027 

3.6 

18.598 

686 

HP 

0.272 

40.6 

0.748 

18.591 

696 

1.091 

107.5 

1.973 

18.546 

571 

1.091 

108.9 

1.958 

18.543 

483 

2.227 

159.0 

2.814 

18.498 

526 

3.786 

198.0 

3.278 

18.452 

416 

6.235 

247.0 

3.945 

18.407 

329 

8.795 

282.8 

4.917 

18.349 

242 

10.354 

293.0 

4.931 

18.313 

236 

11.022 

235.0 

4.329 

18.399 

185 

LML 

12.526 

239.0 

4.337 

18.374 

167 

0.273 

29.9 

0.601 

18.637 

734 

Hr 

1.135 

93.3 

1.880 

18.604 

653 

2.547 

129.4 

3.256 

18.555 

543 

Mac 

5.761 

173.2 

3.950 

18.483 

479 

9.297 

180.3 

4.071 

18.427 

140 

*Amplifier  Impedance 
HP  - < 1 mtt 

Mac  - for  ring  A 100  fi,  for  ring  B 66  ft 
CML  - 64  n 


Close  scrutiny  of  the  data  for  ring  A in  table  1 shows  that  inter- 
changing the  5-W  Hewlett-Packard  (HP)  amplifier  and  200-W  Macintosh 
(Mac)  amplifier  had  very  little  effect  on  the  steady-state  variables 
(current,  mechanical  displacement,  and  resonance  frequency) . However, 
estimated  values  for  Qm  determined  from  the  current  decay  after  vol- 
tage cutoff  differed  by  about  15  percent.  This  discrepancy  indicates 
that  another  potential  problem  in  determining  Qm  values  from  current 
decay  information  exists  because  of  the  influence  of  the  amplifier  im- 
pedance, which  becomes  part  of  the  decay  circuit  during  the  off-time 
of  the  driving  pulse.  This  particular  type  of  problem  may  be  circum- 
vented by  using  the  previously  mentioned  substitution  method  --  namely, 
inserting  an  appropriate  resistance  for  the  ceramic  ring  in  order  to 
obtain  an  equivalent  input  current  value  at  the  same  driving  (reso- 
nance) frequency  and  drive  level,  and  then  computing  the  variation  in 
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the  quality  factor  in  terms  of  the  low-level  value  (Q^  = M / R«) , 

assuming  the  mass  reactance  term  to  be  essentially  constant  and  the 
blocked  impedance  of  the  ring  to  be  much  larger  than  the  mechanical 
resistance  at  resonance.  Appendix  A gives  a brief  description  of  the 
influence  of  the  amplifier  impedance  on  the  quality  factor  in  a sim- 
plified electrical  circuit  and  demonstrates  the  need  for  considering 
such  a potential  problem  with  current-decay  type  measurements. 

The  apparent  decrease  in  the  current  and  displacement  of  ring  A 
when  the  5-kVA  Communications  Measurement  Laboratory  (CML)  amplifier 
was  used  was  not  due  to  the  amplifier  itself  but,  rather,  to  a thin 
sheet  of  cellophane  that  was  placed  between  the  lower  electrode  of  the 
ring  near  the  electrical  lead  and  the  grounded  aluminum  housing  of  the 
gas  suspension  system  in  order  to  prevent  arcing  at  the  higher  drive 
levels.  The  cellophane  added  a small  amount  of  stiffness  that  resulted 
in  an  increase  in  the  resonance  frequency.  The  effect  of  the  cello- 
phane on  the  measured  data  was  taken  into  account  by  extending  the 
curves  in  the  region  of  saturation  --  as  demonstrated  in  figure  11, 
which  is  a plot  of  the  steady-state  strain  S^,  determined  from  the 
measured  displacement  ur=b  (where  b is  the  outside  radius  of  the  ring) 
and  input  current  I,  versus  the  electric  field,  all  in  rms  values. 

The  last  points  on  the  graph  (at  E = 14  kV/m)  correspond  approximately 
to  the  field  level  at  which  ring  A fractured.  The  values  for  and 
I at  this  field  level  are  not  included  in  table  1 because  of  the  bi- 
zarre behavior  of  the  strain  value,  which  was  attributed  to  faulty 
positioning  of  the  displacement  probe.  It  should  be  noted,  however, 
that  the  last  three  points  for  the  current  track  well  and  appear  to 
have  a linear  relationship.  Therefore,  the  strain  value  at  fracture 
is  estimated  from  the  extended  strain  curve  at  the  critical  field  to 
have  a value  = 159  x 10  6 m/m.  This  value  was  used  to  determine 
the  stress  level  at  which  fracture  occurred. 

LINEAR  INTERPRETATION 

Inasmuch  as  there  is  presently  no  precise  nonlinear  mathematical 
theory  for  predicting  the  high-strain  behavior  of  ceramics  at  sonic 
frequencies,  the  experimental  data  presented  in  this  report  are  inter- 
preted solely  in  terms  of  the  variation  of  the  low-level,  linear  ce- 
ramic parameters  as  a function  of  state  variables  (electric  field, 
stress,  etc.).  The  definitions  given  to  the  low-level  parameters  (or 
combinations  of  parameters)  are  determined  from  linear  theory,  starting 
with  the  linear  equations  of  state  for  a mechanical  one-dimensional 
vibrator.  For  the  rings  under  consideration  in  this  investigation, 
for  instance, 

S1  * S11T1  * d31E3  md  (2a) 


w. 
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D3  ' d31Tl  * e33b3  (2I» 

relate  strain  Sj  and  electric  displacement  D3  to  stress  Tj  and  elec- 
tric field  £3,  where 

E 

s . is  the  compliance  coefficient  at  constant  electric  field 
(Ej  = E2  = 0), 

T 

is  the  permittivity  at  constant  stress  (T^  = T^  = 0) , and 
^31*  d31  are  Piezoelectric  coefficients. 


The  coefficients,  which  are  constants  in  the  linear  case,  are  allowed 
to  become  level -dependent  when  the  above  equations  are  applied  to  the 
large-signal  case  (nonrigorous  mathematical  approach).  Reciprocity 
should  not  be  assumed  when  nonlinear  effects  are  present;  therefore, 
the  prime  on  d3j  in  equation  (2b)  indicates  that  it  may  have  a value 
different  from  d3j  in  equation  (2a) . 


The  dynamic  equations  for  the  ring  are  obtained  by  combining  the 
equations  of  state  with  the  expression  for  the  stress  due  to  inertia 


-r  2 2C 

T = pw  r S, 
lml 


(3) 


where 


r^  is  the  mean  radius  of  the  ring, 
p is  its  density,  and 
u is  the  angular  frequency. 

an  expression  for 

(4) 


Substituting  equation  (3)  into  equation  (2a)  yields 
the  dynamic  strain  -- 


S1  = 


d31E3 


r,  22  ei 
L1  • “ rm  psllJ 
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In  reality,  the  compliance  coefficient  is  complex, 

E E'  . E" 

S11  " S11  " JSH  ’ 

and  its  components  are  utilized  in  the  expressions  for  the  resonance 
frequency 

E / 2 E’  \’1/2  f5-, 

“r  = PS11  j * (5) 

and  the  mechanical  quality  factor  at  resonance, 

. - E'  / E" 

Si  ' S11  / S11 

In  terms  of  these  parameters  the  dynamic  strain  becomes 


si  = S = 
1 m 


d31E3 


Sm  = "jd31QME3  at  “ = “x 


where  Sm  is  the  mean  strain  in  the  ring  defined  as  Sm  = ur=r  jr 

The  input  current  to  the  ceramic  ring  may  also  be  determined  by 
manipulation  of  equation  (2a)  into  equation  (2b)  such  that 

°3  = j^A  = d31  ( 4"  S1  ' ~~T~  E3)  + E33§  ({ 

\ S11  S11  / 


since  q - /D3dA  and  I = dq/dt  = juq,  where  q is  the  total  electrical 
charge  and  A is  the  electrode  area  of  the  ceramic  ring.  Rearrangement 
of  equation  (8)  in  conjunction  with  equation  (6)  yields 
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I 

jwA 


d31d31 


’ll 


. T /,  d31d31  \ 

C33  I1  T E ) 
y e33Sll  j 


(9) 


At  resonance,  w = io£,  equation  (9)  reduces  to 


I = 


u AQ,,d71  di. 
r XM  31  31 


E3« 


’ll 


Cio) 


which  is  simply  the  motional  component  of  the  total  input  current 
under  the  assumption  that  the  blocked  component  of  the  current  is  neg- 
ligible (that  is, 

E.  T 
u Ae_, 
r 33 

see  appendix  B) . 


It  is  apparent,  then,  that  figure  11  is  simply  a plot  of  the  mag- 
nitudes of  equations  (7)  and  (10)  at  resonance  of  the  ring,  with  the 
mean  strain  values  replaced  by  the  strain  values  at  the  outside  radius 
of  the  ring  because  these  latter  values  were  more  directly  related  to 
the  actual  measurements.  The  relationship  between  the  strain  at  r = b 

and  the  mean  strain  at  r = r is  given12  as 

m 


r=r 


r=b 


1.153, 


(11) 


E 

where  o is  Poisson's  ratio  for  the  ceramic  material  at  constant 
electric  field  and  has  been  given  the  value  of  0.33  for  this  exercise. 
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The  linear  low-level  piezoelectric  coefficient  d3i  was  determined 
from  the  graph  (based  on  a value  of  = 696  (E  = 0.272  kV/m))  to  be 
137  x lO-12  m/V,  which  compares  well  with  the  value  130  x 10 ~*2  m/V 
calculated  from  the  low-level  admittance  plot  of  ring  A.  Figure  12  is 
a similar  type  of  plot  for  ring  B;  the  linear  piezoelectric  coefficient 
for  this  case  was  determined  to  be  d3^  = 118  x 10"12  m/V,  compared  with 
the  value  of  109  x 10" 12  m/V  from  the  admittance  plot  of  ring  B at  a 
similar  drive  level. 

Figure  13  is  a composite  of  the  strain  versus  electric  field 
curves  from  figures  11  and  12;  it  shows  the  relative  difference  in 
strain  at  fracture  for  the  two  rings,  as  well  as  the  nonlinear  varia- 
tion of  the  product  d3i  Qm,  and  demonstrates  the  need  for  a precise 
measurement  of  Qm  in  order  to  determine  the  behavior  of  d3i  with  drive 
level  at  resonance. 

CALCULATED  PARAMETERS 

Under  the  assumption  that  the  values  of  in  table  1 are  valid, 
equation  (7)  can  be  used  to  evaluate  d3i  at  resonance,  where 


with  the  results  shown  in  figure  14.  The  curves  were  not  drawn  point 
to  point  but  are  the  best  fit  to  the  available  data  points  (haphazard 
at  the  higher  drive  levels,  especially  at  fracture)  and  might  not  rep- 
resent the  true  behavior  of  d3i  with  drive  level.  Nevertheless,  there 
appears  to  be  a definite  tendency  for  the  piezoelectric  coefficient  to 
decrease  with  the  increased  drive  level,  as  shown  by  the  departure  of 
the  data  points  for  ring  A from  the  linear  low-level  value.  The  same 
comment  holds  for  ring  B.  This  result  is  in  disagreement  with  previ- 
ous observations, 13  which  indicated  that  d3i  increased  with  ac  drive 
level  for  similar  types  of  ceramic.  It  is  interesting  to  note,  how- 
ever, that  the  low-level  dynamic  value  of  d3i  was  found  to  decrease14 
with  increased  static  bias  stress,  although  the  magnitude  of  the  de- 
crease was  relatively  small  at  the  comparable  stress  levels  encoun- 
tered here.  It  may  be  that  the  interpretation  given  to  the  high-level 
dynamic  value  of  d3i  should  be  altered  to  account  for  any  losses15 
associated  with  the  piezoelectric  coefficient  at  high  stress  levels. 

If  this  is  done,  then  phase  information  must  be  included  in  the  meas- 
urement program.  Nevertheless,  the  data  presented  by  Gerson10  were 
based  on  the  assumption  that  d3i  remained  constant  with  stress  level 
(this  particular  point  will  be  discussed  in  more  detail). 

Clarification  of  these  discrepancies  can  only  be  achieved  by  more 
detailed  measurements  on  additional  ceramic  specimens. 


22 


(ui/uirfl  •'s  NlVdlS 


Figure  13.  Comparison  of  Strain  versus  Electric  Field  Curves  for  Rings  A and  B 
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Values  for  the  mean  strain  SB  and  Young's  modulus  1/sji  were 
calculated  for  each  electric  field  point  in  table  1 by  using  equations 
(11)  and  (5),  and  the  resulting  data  are  listed  in  table  2 along  with 
the  mean  stress 

T = p<u2r2S  = — 4t-  S at  u = df  . (13) 

m m m E'  m r 

S11 

The  assumed  densities  were  p = 7.54  x 10 3 kg/m3  for  ring  A and  7.65  x 
10 3 kg/m3  foT  ring  B. 


Table  2.  Calculated  Values  for  Ceramic  Rings 


Mean 

Strain 

Sm 

(um/m) 


Electric 

Young ' s 

Ampli- 

fier* 

Field 

e3 

Modulus 

F 1 

1/sfi 

(kV/m) 

(1010N/m2) 

0.027 

8.407 

HP 

0.272 

8.401 

1.091 

8.360 

1.091 

8.357 

2.227 

8.317 

Mac 

3.786 

8.275 

6.235 

8.235 

8.795 

8.183 

10.354 

8.151 

CML 

11.022 

8.228 

12.526 

8.206 

HP 

0.273 

8.565 

1.135 

8.535 

2.547 

8.490 

Mac 

5.761 

8.424 

9.297 

8.373 

71.13 

102.21 

119.08 

143.30 

178.62 

179.14 

157.26 

157.54 


21.83 

68.30 

118.28 

143.48 

147.89 


Mean 

Stress 

r 

m 

(lbf/in? 
x 10 3) 

(MN/m2) 

0.33 

2.28 

0.87 

6.00 

0.86 

5.93 

1.23 

8.48 

1.43 

9.86 

1.71 

11.79 

2.12 

14.62 

2.12 

14.62 

1.88 

12.96 

1.87 

12.89 

0.27 

0.85 

1.86 

5.86 

1.45 

10.00 

1.76 

12.14 

1.79 

12.34 

lifier  Impedance 


- < 1 m£l 

Mac  - for  ring  A 100  fi,  for  ring  B 66  fl 
CML  - 64  fi 


TR  5263 


The  mechanical  quality  factor  Q^,  the  proportional  decrease  in 
Young's  modulus  (Y^i  = l/sf{),  and  the  proportional  decrease  in  the 
piezoelectric  strain  coefficient  d3i  are  plotted  versus  the  mean  stress 
(rms)  values  in  figures  15  and  16  for  rings  A and  B.  Here  Qm  and  l/s^i 
follow  the  same  general  trend  of  variation  as  Gerson's  measurements, 
and  the  magnitudes  of  the  decreases  for  these  quantities  are  strikingly 
similar  to  his  results  if  aging  of  the  rings  is  allowed  for  and  the 
stress  values  quoted  in  this  report  are  converted  to  peak  values. 
Gerson's  observations  that  the  values  increase  with  aging  and  that 
the  variation  with  stress  level  becomes  more  linear*  are  confirmed  by 
the  NUSC  measurements.  Gerson  also  found  that  the  variation  in  Young's 
modulus  with  drive  level  was  less  severe  for  well -aged  ceramic  speci- 
mens, and  this  is  confirmed  by  the  measurement  data  presented  in  this 
report.  Minor  differences  in  values,  as  measured  by  Gerson  or  at  NUSC, 
can  be  attributed  to  differences  in  the  ceramic  sample  properties, 
which  are  outlined  in  table  3 for  comparison. 


Table  3.  Properties  of  Ceramic  Samples 


Ceramic 

Property 

Gerson1 °t 

NUSC* 

Ring  A 

Ring  B 

k31 

0.247 

0.355 

0.326 

K^dkHz) 

1240 

1278 

1077 

d31 

-92xl0‘12  m/V 

-130xl0'12  m/V 

-109xio'12  m/V 

ye 

11 

7.92xl010  N/m2 

8.39xl010  N/m2 

8. 57x10 10  N/m2 

Qm 

530 

833 

965 

tPbo.91*  Sro#o6  Zr0.53  Tio.47  O3  ceramic,  radially  poled  rings,  sus- 
pended by  thread, *and  measured  200  days  after  poling. 
tNavy  Type-I  ceramic,  axially  poled  rings,  suspended  by  rubber 
band,  and  measured  3 years  after  poling. 


•The  plots  of  Qm  versus  peak  dynamic  stress  in  Physical  Acoustics1 
(p.  217)  are  apparently  for  newly  poled  ceramics  inasmuch  as  the 
initial  slope  of  the  curves  is  quite  steep. 
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Figure  15.  Mechanical  Quality  Factor,  Decrease  in  Young's  Modulus, 
and  Decrease  in  Piezoelectric  Coefficient  versus  Mean  Stress  for  Ring  A 
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Figure  16.  Mechanical  Quality  Factor,  Decrease  in  Young's  Modulus, 
and  Decrease  in  Piezoelectric  Coefficient  versus  Mean  Stress  for  Ring  B 
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The  variation  of  (I31  with  stress  is  quite  severe  in  figures  15  and 
16,  and  there  appears  to  be  a region  where  the  change  becomes  maximum 
--  although  this  observation  cannot  be  considered  factual  inasmuch  as 
the  region  is  one  of  uncertainty  (very  near  fracture).  Nevertheless, 
the  information  definitely  indicates  that  d3i  decreases  at  relatively 
low  drive  levels  (refer  also  tq  figure  14) , and  this  result  appears  to 
conflict  with  Gerson's  assumption  that  d3i  remains  constant  with  drive 
level.  It  should  be  remembered,  however,  that  the  equation  used  by 
Gerson  to  determine  the  stress  from  the  measured  current  is  the  same  as 
equation  (2b),  if  the  free  component  of  the  current  (aipAc^Es)  is  ig- 
nored, where  d3j  has  been  replaced  by  d^i  to  account  for  nonlinear 
effects.  The  stress  values  used  in  the  present  observations  were  ob- 
tained from  directly  measured  quantities  (equation  (13)),  and  the  re- 
sulting stress  values  were  used  in  conjunction  with  the  measured  cur- 
rent values  to  compute  the  behavior  of  the  piezoelectric  stress 
coefficient  d$i  with  stress  level  as  depicted  in  figure  17  (both  ab- 
scissa and  ordinate  are  in  rms  values)  . This  plot  indicates  that  d^i 
behaves  quite  well  with  drive  level,  with  slight  variations  in  value  at 
high  stress  level,  and  that  the  initial  low  level  value  is  in  very  good 
agreement  with  the  piezoelectric  strain  coefficient  d3i  determined  from 
the  strain  versus  electric  field  relationship  in  figures  11  and  12. 
Therefore,  the  results  presented  by  Gerson  are  still  valid  as  long  as 
it  is  understood  that  it  is  d^  (C/N)  that  is  essentially  constant  with 
stress  level  and  not  d3i  (m/V) , whose  values  are  equal  only  for  the 
low-level  linear  case. 

It  is  interesting  to  note  that  the  variation  in  d^i  (low  current 
and  high  stress)  for  ring  A follows  the  same  general  trend  as  the  pie- 
zoelectric strain  coefficient  d3j  measured  by  Woollett16  under  snail 
strain  and  high  field  conditions  — determined  from  equation  (2a)  by 
driving  long,  thin,  thickness-poled  bars  at  a frequency  well  below  the 
fundamental  length  resonance  of  the  bars  where  the  stress  Tj  may  be 
neglected.  The  implication  here  is  that  different  kinds  of  driving 
conditions  for  the  ceramic  might  necessitate  the  use  of  different  co- 
efficients in  the  equations  of  state  if  these  equations  are  used  to  ex- 
press the  behavior  of  ceramic  materials  under  nonlinear  conditions. 

For  example,  d^j  (low  current  and  high  stress  in  C/N)  might  be  equiva- 
lent to  d3i  (low  strain  and  high  electric  field  in  m/V) , whereas  d3i 
(low  electric  field  and  high  strain)  might  be  equivalent  to  d^i  (low 
stress  and  high  current).  If  such  a situation  exists,  then  different 
symbols  would  probably  be  needed  for  the  different  sets  of  conditions.* 


*No  proof  exists  at  present  to  substantiate  such  comments,  which  are 
made  solely  to  activate  the  imagination  of  researchers  in  transducer 
ceramic  studies. 
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Figure  18  is  a plot  of  the  magnitude  of  equation  (10)  with  sfj  re- 
placed by  the  real  part  sff  of  its  complex  quantity  (i.e.,  1/Qm  was 
considered  to  be  much  less  than  unity  for  all  practical  considerations) 
although  it  does  not  offer  any  additional  information,  it  is  included 
for  completeness. 


Figure  18.  Behavior  of  Product  Q^^d^  for  Rings  A and  B 

The  maximum  breaking  (fracture)  strength  of  the  specimens  tested 
was  computed  by  using  Sinsky's12  correction  for  determining  the  max- 
imum strain  at  the  inside  radius  of  the  rings  from  the  strain  at  the 
outside  radius  of  the  rings  (the  correction  for  the  dimensions  of  the 
rings  used  here  and  for  a*  = 0.33  was  Sa  = (1.345)51,)  and  the  value  of 
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Young's  modulus  at  fracture.  Calculations  for  ring  B were  straightfor- 
ward, but  those  for  ring  A had  to  be  corrected  for  the  influence  of  the 
cellophane  insulator  at  the  higher  drive  levels.  Figure  19  is  a plot 
of  the  resonance  frequency  of  the  rings  versus  the  electric  field  level 
and  was  used  to  estimate  the  true  resonance  frequency  at  fracture  for 
ring  A (f|  = 18.255  kHz),  from  which  Young’s  modulus  was  calculated  by 
using  equation  (5).  The  fracture  strain  (159  pm/m)  for  ring  A was 
taken  from  the  extended  portion  of  the  strain  curve  in  figure  11,  and 
that  for  ring  B (S^  = 128  pm/m)  was  taken  from  the  strain  curve  in  fig- 
ure 12.  The  maximum  (peak)  dynamic  stress*  was  computed  by  using  equa- 
tion (13),  with  the  mean  values  replaced  by  the  values  at  r=a  to  yield 
3550  lbf/in. 2 (24.5  MN/m2)  and  2950  lbf/in.2  (20.3  MN/m2)  for  rings  A 
and  B,  respectively.  If  the  low-level  Young's  modulus  is  used  in  the 
computations,  the  values  are  slightly  higher  --  3690  lbf/in.2  (25.4 
MN/m2)  and  3030  lbf/in.2  (20.9  MN/m2).  The  computed  stress  limits  are 
in  agreement  with  Gerson's  observations  and  with  the  dynamic  values 
quoted  by  commercial  manufacturers, 17  and  are  only  10  to  15  percent 
lower  than  the  average  dynamic  fracture  strengths  quoted  by  Smith  and 
Rice18  for  similar  types  of  ceramic  material  (for  Edo  Western  4100 
lbf/in.2  (28.3  MN/m2)  and  for  Channel  3300  lbf/in.2  (22.8  MN/m2)). 

Other  observers13  have  claimed  values  as  high  as  5000  to  6000  lbf/in.2 
(34.5  to  41.4  MN/m2)  for  the  peak  dynamic  tensile  strength  of  similar 
hard  lead-zirconate-titanate  ceramic  supplied  by  a different  manufac- 
turer. All  of  the  values  quoted  might  conceivably  be  correct  inasmuch 
as  the  ultimate  stress  level  depends  on  such  variables  as  manufacturing 
processes,  age  after  polarization  when  the  ceramic  is  tested,  degree  of 
pore  concentrations  in  the  ceramic,  and  heating  effects  at  high  drive 
level.  Clarification  of  the  influence  of  these  variables  on  fracture 
strength  can  only  be  achieved  through  further  measurements. 

It  should  be  noted  that  in  all  cases  the  dynamic  fracture  strength 
of  Navy  Type- I ceramic  is  about  50  or  60  percent  of  the  static 
strength.13 


*It  should  be  noted  that  the  rings  actually  broke  at  the  next  higher 
field  level  tried  in  the  measurements  but  that  the  last  observed  field 
data  point  was  the  one  used  for  stress  fracture  calculations.  The 
quoted  fracture  levels  are  still  valid,  however,  since  fracture  occurs 
in  a region  of  strain  (or  stress)  saturation. 
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Figure  19.  Resonance  Frequency  versus  Electric  Field 
for  Rings  A and  B 

NONLINEAR  THEORETICAL  ANALYSIS 

The  lack  of  an  existing  precise,  practical,  nonlinear  mathmetical 
theory  for  predicting  the  behavior  of  piezoelectric  ceramics  in  sonar 
transducers  has  been  noted.  The  key  word  is  sonar  (sonic,  as  opposed 
to  ultrasonic,  frequencies),  since  ceramic  properties  useful  to  trans- 
ducer design  engineers  are  usually  determined  on  specific  geometric 
specimens  (see  figure  20)  driven  electrically  at  their  fundamental  res- 
onance frequencies  and  sometimes  at  frequencies  well  below  the  funda- 
mental resonance. 

One  phase  of  this  investigation  was  to  develop  a simplified  non- 
linear theory  that  could  be  used  to  interpret  the  experimental  data 
collected.  A simplified  mathematical  approach  is  currently  preferred 
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as  a tool  because  of  the  complexity  that  would  be  involved  if  a 
rigorous  continuum-physics  approach5  were  used  to  analyze  composite 
transducer  structures  and  transducer  arrays.  Even  a simplified  ap- 
proach to  nonlinear  theory  involves  considerable  complication,  but  with 
appropriate  assumptions  and  boundary  conditions  usable  equations  can  be 
developed. 


ELECTRODES  ON 


2b  » h > (b  • a)  I » w > h I » w > h 


Figure  20.  Typical  Ceramic  Test  Specimens 


APPROACH 


An  approach  similar  to  that  used  for  determining  third-order  pie- 
zoelectric coefficients20  from  ultrasonic  measurement  information  will 
be  followed.  The  same  approach  can  not  be  used  because  at  ultrasonic 
frequencies  equations  of  state  are  usually  written  in  terms  of  the  elas- 
tic stiffness  coefficients2^  c^  at  constant  electric  field  or  c°  at 
constant  electric  displacement  (which  are  the  appropriate  parameters 
measured  with  ultrasonic  techniques),  whereas  the  equations  of  state 
applicable  to  transducer  theory  (sonic  frequencies)  are  usually  written 
in  terms  of  the  elastic  compliance  coefficients  s^  or  s°.  Although 
there  exists  a relationship  between  s and  c (with  the  appropriate 
superscripts) ,22 


s. . = 


[-1)1+J  dcij 
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where  Ac  is  the  determinant  of  the  Cjj  terms  and  Acjj  is  the 
minor  obtained  by  suppresssing  the  i-th  row  and  j-th  column,  there  ap- 
pears to  be  some  doubt  as  to  its  rigorousness  if  used  to  determine  low 
frequency  coefficients  from  the  measured  high  frequency  coefficients. 

In  addition,  equation  (14)  was  derived  from  the  linear  constitutive 
equations  of  state  and  the  problem  being  addressed  here  is  a nonlinear 
one. 

Therefore,  the  expanded  nonlinear  constitutive  equations  of  state 
for  a piezoelectric  medium  (keeping  terms  through  cubic  in  the  poten- 
tial energy  density)  and  the  field  equation  of  motion  for  an  elastic 
continuum  were  applied23  to  a thin-walled  ceramic  ring  (see  figure  20) 
operating  near  its  fundamental  hoop-mode  resonance.  This  specific 
cylindrical  geometry  was  chosen  for  ease  of  applying  nonlinear  theory 
as  well  as  for  the  experimental  aspect  of  this  investigation,  which  in- 
volved measurements  of  ceramic  rings.  The  assumptions  associated  with 
a one-dimensional  ring  vibrator  are  that  all  stress  components  are 
identically  zero  except  the  tangential  component  T^,  and  all  variables 
are  independent  of  the  axial  and  tangential  coordinates  z(3)  (the  di- 
rection of  polarization  in  the  ceramic  ring)  and  0(1),  respectively; 
with  an  abbreviated  form  of  tensor  notation21  the  reduced  equations  for 
the  strain  m and  the  electric  displacement  D3  become 


nl  = S11T1  + d31E3  + 2 S111(T1)  + 2 d331(E3) 


D3  = d31Tl  + E33E3  + 2 d311(1V  + 2 e333(E3)  + d331E3Tl’ 


where 

E T 

second-order  coefficients  sn,  d3i,  and  £33  are  the  usual  linear 
elastic,  piezoelectric,  and  dielectric  coefficients, 

E T 

the  remaining  coefficients  sm,  d33i,  d3n,  and  £333  are  the  non- 
linear or  third-order  coefficients  of  the  piezoelectric  ceramic, 
and 


is  the  (17) 


j 


approximate  equation  of  motion  for  the  ceramic  ring. 
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The  Lagrangian  strain  component  m 


in  the  tangential  direction  is  given 


(18) 


When  the  nonlinear  terms  are  ignored,  equations  (15)  through  (18)  re- 
duce to  the  common  linear  expressions.  Reciprocity  should  not  be  as- 
sumed when  nonlinear  effects  are  present,  but  not  doing  so  would  only 
aggravate  the  already  complex  situation.  The  first  step,  therefore,  is 
to  see  if  a useful  wave  equation  can  be  derived  from  the  four  nonlinear 
equations  for  a thin-walled  ring. 


When  the  ring  is  excited  very  near  resonance,  large  stresses  and 
strains  are  induced  in  the  ring  by  low  electric  field  levels  because  of 
the  high  of  the  material.  The  controlling  parameters  for  this  case 
are  the  elastic  coefficients,  and  only  the  linear  electric  field  term 
will  be  retained.*  Under  these  conditions  the  stress  becomes  approx- 
imately23 


V-r  1 +J1  * 2 

sm  I y 


nl  ' d31E3| 


(19) 


With  the  use  of  equation  (18)  the  radical  can  be  expanded  (keeping 
terms  to  the  order  consistent  with  previous  approximations)  to  yield 

2 


T1  “ 


1 

E 

S11 


u 


1 - 


2s 


11 


’ll 


31 


’ll 


(20) 


The  last  expression,  in  conjunction  with  equation  (17),  yields  the  non- 
linear wave  equation 


7V  * (»r)  2(1  - ■ (»r)  2 d31E3  - <21> 


♦Note  that  equations  (15)  and  (16)  may  also  be  used  at  very  low  fre- 
quencies (Tj«0)  to  obtain  information  about  the  third-order  piezoelec 
trie  (d33i)  and  dielectric  (C$33)  coefficients. 
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where 


n = rr-  ■ , and 
m 


r i rm  (the  mean  radius  of  the  ring) . 


The  solution  of  equation  (21)  is  by  no  means  straightforward,  if  even 
possible.  Two  different  approaches  were  tried  in  order  to  obtain  in- 
formation from  the  nonlinear  wave  equation.  First,  an  approximate  per- 
turbation for  the  strain  variable  n was  tried,  and,  second,  an  analog 
computer  approach  was  used  to  obtain  a solution. 

PERTURBATION  SOLUTION 


It  was  initially  assumed  that  the  strain  variable  n could  be 
written  as  the  sum  of  two  parts 


n 


,<°>  ♦ n(1) 


(22) 


where 


is  the  solution  of  the  linear,  unperturbed  (y  = 0)  equation 
and 

is  the  remaining  or  perturbation  part  n^/n^  |<<1^. 

If  it. is  assumed  also  that  the  electric  field  is  sinusoidal  (E3 
= fc'3ejft)t),  then  q(0)  and  qU)  can  be  written  as 


n 


(0) 


iW*  and  n(1) 


m 


nCl)ej2“t 

m 


(23) 
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where  the  maximum  value  of  the  fundamental  component  is 

(0)  _ d31g3 

\ ~ r ■» 


and  the  maximum  value  of  the  second  harmonic  is 


(24) 


(25) 


In  order  to  keep  the  strain  components  finite  it  is  necessary  to  con- 
sider the  linear  compliance  coefficient  sfi  as  a complex  quantity. 
Then  equations  (24)  and  (25)  become,  at  w = 


= “3^31^63  (see  equation  (7))  and 


(26) 


assuming  that  Qm  is  large  (4/Qm«3).  Therefore,  it  would  appear  that 
spectral  analysis  measurements  of  the  relative  level  of  second  harmonic 
and  absolute  fundamental  displacement  components  would  yield  informa- 
tion about  the  third-order  elastic  compliance  coefficient  sf u . 


In  addition,  the  current  at  resonance  can  be  derived,  retaining 
the  second-order  electric  field  term  in  equation  (16),  as 


I = 


2 E. 
d3l“rA 


’ll 


Qm  + j“rAe 


2“rAd31l4 


ft)' 


* i“rAtl33  1 E3  • 


(27) 
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where  6333  has  been  defined  as 

,2 


s T 
333  ' e333 


1 + 


d311d31 


’333 


& 


d31^d311 


'333 


(•y 


2d 


1 + 


31  11 


311 


H>) 


(28) 


Note  that  the  leading  (linear)  term  in  equation  (27)  is  the  same  as 

4. 


equation  (9)  when  d^  = d31  and  w = co^,  since 


'33 


= £jj{1  - [d,,  / (e 


31 


T E 
'33S11 


)]>. 


Since  it  is  shown  in  appendix  B that  the  fundamental  blocked 
component  (1^,  = ju^AefaEs)  of  the  input  current  is  negligible,  for  all 
practical  purposes,  for  ring  measurements  described  in  this  report,  it 
will  be  assumed  that  its  counterpart  in  the  second-harmonic  component 
is  also  negligible.  Under  these  conditions  the  current  becomes 


I 


d32l"fon 

4 


2“rAd311d31%.  2 

E3 


(29) 


which  demonstrates  that  spectral  analysis  could  be  used  to  yield  infor- 
mation about  the  third-order  piezoelectric  coefficient  inasmuch  as 
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where  l(l)  and  I (0)  are  the  magnitudes  of  the  second-harmonic  and  fun- 
damental current  components,  respectively. 


i-pectral  analysis  measurement  of  the  radial  displacements  and  in- 
put driving  currents  for  the  rings  is  difficult  with  conventional 
equipment  because  heating  problems  force  the  rings  to  be  excited  with 
low-duty- cycle  pulses.  Reliable  measurement  data  would  probably  re- 
quire sophisticated  techniques,  such  as  fast  Fourier  transform  (FFT) 
processing.  This  aspect  has  not  been  thoroughly  considered  because  it 
was  felt  that  perturbation  theory,  although  art  improvement  over  linear 
theory  and  useful  in  explaining  some  aspects  of  nonlinear  behavior, 
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... 


would  not  be  sufficient  to  explain  the  gross  nonlinear  effects  experi- 
enced in  the  ring  measurements.  Therefore,  the  decision  was  made  to 
try  an  alternative  approach  to  the  problem. 

ANALOG  COMPUTER  SOLUTION 


A preliminary  analog  computer  solution  was  sought24  in  which  the 
resulting  response  pulse  shapes  and  jump  phenomena  would  duplicate  ac- 
tual experimental  observations.  Duplication  of  complete  ring  behavior 
was  not  attempted  because  of  insufficient  time.  Instead,  the  analog 
computer  was  programmed  by  using  experimentally  determined  parameters 
from  measurements  of  a single  ring  at  a single  driving  field  in  order 
to  demonstrate  the  capabilities  of  the  analog  computer. 

The  first  predicament  encountered  was  the  necessity  to  incorporate 
some  appropriate  loss  term  in  the  wave  equation  to  keep  n finite. 

Since  the  loss  parameter  in  a simple  mechanical  system  is  associated 
with  the  velocity  component,  it  has  been  assumed  that  a term  propor- 
tional to  the  time  derivative  of  the  second  term  on  the  left-hand  side 
of  equation  (21)  would  account  for  nonlinear  losses.  Hence,  equation 
(21)  becomes 

2 

+ 6(1  - 2yn)  g£  + (<^)2(1  - Yn)n  = (Jb2d31E3  . (31) 


where 


E 


is  the 
linear 


and 


low-level  mechanical  quality  factor  measured  under 
conditions. 


When  the  nonlinear  elastic  parameter  y is  zero,  equation  (31)  reduces 
to  the  conventional  form  found  in  textbooks  on  mechanical  vibrations. 

Equation  (31)  was  simulated  on  an  Electronics  Associates,  Incor- 
porated (EAI)  680  analog  computer  with  the  results  shown  in  figure  21. 
Figure  21a  shows  measured  pulse  shapes  for  an  actual  ring;  the  upper 
trace  is  the  applied  voltage  pulse,  the  middle  trace  is  the  resulting 
current  pulse,  and  the  bottom  trace  is  the  resulting  radial  displace- 
ment of  the  ring.  The  computer  simulation  at  the  same  drive  level  is 
shown  in  figure  21b,  where  the  upper  trace  is  the  driving  pulse  and  the 
lower  trace  is  the  strain  response  to  the  driving  pulse.  Note  that  the 
displacement  pulse  in  figure  21a  has  a strange  shape.  In  all  other 
measurements  made  on  ceramic  rings  during  this  project  it  was  found 
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that  the  displacement  pulse  shapes  were  almost  identical  to  the  current 
pulse  shapes  at  the  same  drive  levels.  Hence,  it  seems  reasonable  to 
conclude  that  there  was  an  instrumentation  problem  with  the  displace- 
ment probe  in  this  particular  case  and,  therefore,  that  comparison 
should  be  made  between  the  analog  strain  pulse  shape  and  the  experimen- 
tal current  pulse  shape.  Further  clarification  of  this  point  could  not 
be  made  because  the  ring  used  for  this  study  had  fractured. 


Figure  21a.  Actual  Figure  21b.  Computer 
Voltage,  Current,  Generated  Voltage  and 
and  Displacement  Strain- Response 


Figure  21.  Comparison  of  Computer  - Generated 
Data  with  Experimental  Data 


It  is  interesting  to  note  that  all  of  the  response  pulses  had  rip- 
ples (similar  to  beat  phenomena)  that  could  not  be  eliminated  by  vary- 
ing the  frequency  and  that  appeared  to  decrease  exponentially  with  time. 
This  phenomenon  was  not  observed  at  the  very  low  driving  fields  (linear 
behavior,  see  figure  7a)  but  became  pronounced  as  the  drive  level  in- 
creased (increased  current  and  strain) . This  same  phenomenon  has  been 
observed  and  discussed  elsewhere25  in  the  field  of  nonlinear  ceramic 
studies . 

Figure  22  depicts  the  real  time  spectrum  of  the  steady-state  por- 
tion of  the  analog  strain  pulse  and  shows  the  presence  of  second-  and 
third-harmonic  components  at  relative  levels  of  -32  dB  and  -43  dB,  re- 
spectively, below  the  fundamental  component.  Information  of  this  na- 
ture could  be  useful  in  estimating  the  relative  magnitude  of  the  second 
harmonic  to  look  for  experimentally  (if  the  results  of  the  perturbation 
theory  presented  previously  are  applied  at  discrete  drive  levels  for 
the  rings) . 

Figure  23  illustrates  the  jump  phenomenon  that  is  so  characteris- 
tic of  nonlinear  problems  and  that  can  be  verified  by  referring  to  the 
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experimental  results  in  figures  7e  through  7g.  Moving  from  left  to 
right  in  figure  23,  the  simulated  strain  (« displacement)  responses  for 
frequencies  slightly  below,  at,  and  slightly  above  resonance  are  depic- 
ted. This  sequence  of  events  demonstrates  a highly  nonlinear  situation 
since  the  magnitude  of  the  strain  is  highly  nonsymmetric  about  reso- 
nance. In  addition,  it  should  be  noticed  that  the  so-called  beat  phe- 
nomenon discussed  previously  is  quite  apparent  here  and  that  it  also  is 
highly  nonsymmetric  about  resonance  — the  effective  number  of  beats 
being  considerably  less  below  resonance  than  at  (and  above)  resonance. 


Figure  22.  Spectrum  of  Simulated  Steady-State 
Radial  Displacement  of  Ring 
(resonance  frequency  f|  = 18.290  kHz) 

It  should  be  reemphasized  at  this  point  that  the  results  presented 
in  the  latter  part  of  this  theoretical  investigation  are  only  prelimi- 
nary in  the  sense  that  the  analog  computer  solution  was  the  last  ap- 
proach attempted  and  there  was  not  sufficient  time  to  completely  evaluate 


i 
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that  approach  to  determine  the  nonlinear  ceramic  behavior.  However,  the 
demonstrated  capability  of  an  analog  computer  simulation  to  reproduce 
actual  experimental  data,  including  highly  nonlinear  behavior  such  as 
beat  and  jump  phenomena,  makes  the  use  of  the  analog  computer,  in  con- 
junction with  derived  theoretical  results  and  experimental  data,  very 
encouraging  for  obtaining  quantitative  information  about  the  nonlinear 
characteristics  of  piezoelectric  ceramics. 

DRIVING  PULSES 


STRAIN  RESPONSES 


f.  = 18,252  Hz  f = 18,290  Hz 

^ o 


f = 18,328  Hz 
u 


Figure  23.  "Junp"  Phenomenon  (simulated) 


CONCLUSIONS  AND  RECOMMENDATIONS 

A measurement  program  was  established  to  probe  the  nonlinear  behav- 
ior of  piezoelectric  ceramic  rings  at  high  strain  levels.  Although  the 
number  of  rings  tested  was  limited,  the  information  obtained  by  using 
linear  piezoelectric  theory  to  interpret  the  experimental  data  (inasmuch 
as  it  was  the  only  analytical  tool  available  at  the  start  of  the  pro- 
gram) included  several  interesting  features.  Among  them  was  the  ac- 
cepted (or  established)  behavior  of  the  mechanical  quality  factor  and 
of  Young's  modulus  (of  the  ceramic)  versus  stress  level.  In  addition, 
initial  measurements  indicated  that  the  piezoelectric  coefficient  d3i 
commonly  used  in  the  equations  of  state  for  piezoelectric  material  does 
not  obey  reciprocity  when  nonlinear  effects  are  present.  It  was  tenta- 
tively found  that  the  strain  coefficient  d3j  (=  S1/E3)  decreased  with 
stress  level,  whereas  d^  (=  D3/T x)  increased  with  stress  level.  These 
observations  need  to  be  verified  by  additional  measurements  on  ceramic 
specimens  supplied  by  other  manufacturers;  it  was  felt  that  the  rings 
measured  in  this  investigation  exhibited  exceptionally  large  nonlinear 
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effects  that  might  not  be  typical  of  Navy  Type- I ceramics.  Also,  the 
measurement  scheme  should  be  improved  to  account  for  phase  information 
and  possible  spectral-analysis  measurements.  The  potential  problems 
associated  with  the  present  measurement  system  need  to  be  better  under- 
stood and  corrected  where  possible  (for  example,  by  improving  the  gas- 
bearing assembly  to  eliminate  unwanted  losses). 

Theoretical  analysis  of  the  nonlinear  behavior  of  a piezoelectric 
ceramic  ring  demonstrated  the  complications  involved  in  finding  a gen- 
eral solution  for  predicting  high-level  nonlinear  behavior  in  ceramics. 
However,  a reduced  nonlinear  differential  equation  developed  for  the 
ceramic  ring  was  programmed  on  an  analog  computer,  and  the  resulting 
outputs  were  strikingly  similar  to  actual  experimental  data.  These 
kinds  of  results,  compared  with  the  complex  aspects  of  a generalized 
nonlinear  theory,  indicate  that  analog  computer  techniques  may  be  the 
only  convenient  means  of  interpreting  and  understanding  ceramic  non- 
linear behavior.  The  analog  computer  solution  must  be  matched  with  ac- 
tual driving  conditions  and  measurement  data  on  ceramic  rings  before 
quantitative  information  can  be  obtained  about  the  nonlinear  elastic 
and  piezoelectric  parameters. 
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APPENDIX  A 


INFLUENCE  OF  AMPLIFIER  IMPEDANCE  ON 


<*M 


The  piezoelectric  ceramic  ring  vibrator,  being  essentially  a cne- 
dimensional  vibrator  near  the  hoop-mode  resonance  frequency,  can  be 
represented  by  an  equivalent  circuit  consisting  of  a blocked  capaci- 
tance Cfc  in  parallel  with  a series  RCL  branch  (see  figure  A-l).  The 
mechanical  quality  factor  for  such  a circuit  is  usually  expressed  in 
terms  of  the  electrical  boundary  conditions.  More  explicitly,  for 
short-circuit  electrical  conditions  (e(t)  = 0) 


“r  = "r  " -f=k  an<i  * -§T a"d 


Jl  c 

1 y y 


wC  R 
r y y 


(A-l) 


for  open-circuit  electrical  conditions  (1=0) 


U)  = 0)  = 

r r 


I.  SCy 
X y (cb  * V 


and 


- k“ 


V 1 - k2 


(A-2) 


where  k is  the  electromechanical  coupling  coefficient  of  the  circuit 
and  is  defined  as 


♦ C 


(A-3) 


The  problem  encountered  in  this  investigation  was  to  determine  what 
happens  to  the  quality  factor  of  the  circuit  when  the  driving  voltage 
e(t)  is  replaced  by  its  equivalent  voltage  source  and  internal  imped- 
ance. This  type  of  problem  is  encountered  when  the  source  voltage  is 
suddenly  terminated  and  the  current  is  allowed  to  decay  through  all  im- 
pedance elements  in  the  circuit,  including  the  source  impedance  — as 
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when  measuring  ceramic  properties  from  current  decay  phenomenon,  which 
is  the  case  at  hand. 


b 


Figure  A-l.  Electromechanical  Equivalent  Circuit  for 
Ceramic  Ring  Vibrator  with  Amplifier 


The  driving  voltage  e(t)  may  be  represented  by  using  Thevenin's 
theorem  for  any  two-terminal  linear  network,  as  a generator  equal  to 
the  open-circuit  voltage  appearing  between  the  terminals  a and  b in 
series  with  an  equivalent  output  impedance  Za.  This  representation 
will  be  used  to  investigate  the  influence  of  the  impedance  of  the  power 
amplifier  (which  is,  in  reality,  a nonlinear  device)  on  the  quality 
factor  of  an  RCL  circuit;  the  scheme  is  commensurate  with  analyzing 
nonlinear  behavior  of  piezoelectric  ceramic  through  linear  equations. 

Amplifier  impedance  is  generally  considered  to  be  a complex  quan- 
tity, whose  magnitude  varies  with  the  operating  frequency  as  well  as 
with  the  load  impedance.  Inasmuch  as  the  exact  components  of  am- 
plifier impedance  were  not  known,  it  was  assumed  that  the  impedance  was 
a real  number  (resistive  component  Ra)  in  order  to  discover  what  in- 
fluence even  this  simple  condition  would  have  on  the  quality  factor  of 
the  circuit;  two  approaches  were  used.  First,  the  quality  factor  was 
determined  by  simply  considering  the  effects  of  an  additional  resist- 
ance in  a simplified  circuit  under  steady-state  conditions.  Second, 
the  quality  factor  of  a similar  circuit  under  "pulse-driving"  condi- 
tions was  determined  from  the  current  decay  in  the  circuit  during  the 
"off-time"  of  the  electrical  driving  signal. 
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Adjustment  of  the  amplifier  impedance  is  most  easily  obtained  in  terms 
of  the  blocked  impedance  of  the  ring  by  letting 


r * 1 £ — I , for  0 < X < ».  (A-7) 

a | Vb  I “ " 

The  end  points  are  easily  verified:  X 0 implies  Ra  -*■  0 and 

Ca  ■+  » [ 1/  (wrCa)  -*■  0]  such  that  the  only  remaining  circuit  elements  are 

Cy,  Ly,  and  Ry,  which  yield  the  short-circuit  electrical  results  stated 

in  equation  (A-l);  and  X ® implies  Ra  0 and  Ca  Cjj  such  that  the 
remaining  circuit  elements  are  Cb,  Cy,  Ly,  and  Ry,  which  yield  the  open- 
circuit  electrical  results  stated  in  equation  (A-2) . 

The  objective  is  to  determine  what  happens  to  the  resonance  fre- 
quency and  quality  factor  of  the  circuit  for  intermediate  values  of  X 
(or  Ra).  The  answer  is  obtained  by  rearranging  and  evaluating  equations 
(A-5)  and  (A-6)  in  terms  of  the  variable  X as  follows  — 
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Evaluation  of  equation  (A-8)  indicates  that  the  resonance  frequency 
u>r  begins  at  ul  and  asymptotically  approaches  w|/  Vl-k^  = oi£  as  X 
increases  positively  from  zero  to  infinity.  For  the  ceramic  ring  case, 
where  k » 0.3,  the  variation  between  extreme  resonance  points  is  only 
5 percent  and,  therefore,  may  be  considered  constant  for  all  practical 
purposes  and  equated  to  unity  in  evaluating  equation  (A-9) . With  this 
assumption,  and  using  values  that  are  representative  of  the  ceramic 
rings  measured  in  this  investigation  (f|  = 20  kHz,  = 500,  k = 0.3, 

and  Cb  * 10-9  F),  equation  (A-9)  is  plotted  in  figure  A-3  as  a function 
of  the  dimensionless  variable  X.  The  plot  is  shown  to  be  symmetrical 
about  X = 1,  but  this  is  due  to  having  neglected  the  slight  increase  in 
the  ratio  u)r/u)|  as  X increases.  The  dotted  portion  of  the  curve  is  in- 
tended to  account  for  the  slight  5 percent  increase  in  Qm  above  Q|  as  X 
becomes  infinite  (to  agree  with  equation  (A-2)) . The  quality  factor 
drops  to  a minimum  (Qm  = 20)  when  X = 1,  or  when  the  amplifier  imped- 
ance equals  the  impedance  of  the  blocked  branch  of  the  ceramic  ring, 
about  8 kfl  at  The  same  conclusion  can  be  deduced  directly  from 

equation  (A-9)  by  simply  maximizing  the  term  in  the  denominator  con- 
taining X (specifically,  differentiate  X/(l+X2)  with  respect  to  X,  set 
the  result  equal  to  zero,  and  solve  for  X)  . 


The  nonlinear  variation  of  Qm  versus  amplifier  impedance,  even 
if  applied  to  the  cases  under  consideration  in  this  report,  would  not 
be  as  critical  as  the  plot  indicates  inasmuch  as  the  amplifier  imped- 
ance was  much  less  than  the  blocked  impedance  of  the  rings.  This  fact 
was  verified  by  making  some  quick  measurements  with  the  200-W  Macintosh 
amplifier.  The  66-fi  tap  on  the  output  of  the  amplifier  was  used  and 
corresponded  to  the  expected  impedance  of  the  rings  at  resonance.  The 
open-circuit  voltage  eoc(t)  and  the  voltage  e(t)  across  several  load 
resistors  Rl  were  measured  at  three  discrete  frequencies,  and  the  magni- 
tude of  the  amplifier  impedance  |Za|  was  determined  from  the  equation 

Rt 

I e(t)  | = | e (t)  | ^ , where  R h |Z  |.  (A-10) 

(RL  + Ra) 

The  measurements  were  made  by  using  continuous  wave  excitation,  and  the 
magnitudes  of  the  voltages  were  determined  just  before  clipping  appeared 
in  the  output  signal.  The  data  are  listed  in  table  A-l  and  indicate 
that  at  the  ring  operating  frequency  of  18.5  kHz,  which  is  near  the  up- 
per limit  of  the  amplifier,  the  amplifier  impedance  is  approximately 
20  fi  for  a load  resistance  of  100  fi.  If  this  value  were  assumed  to  be 
a real  number,  independent  of  frequency,  as  was  assumed  for  figure  A-3, 
then  the  corresponding  Qm  would  be  445,  or  90  percent  of  the  short-cir- 
cuit value.  The  measurements  also  indicate  that  the  correction  would 
be  even  less  for  the  rings  at  high  drive  levels  because  | Za| decreases 
with  inc -eased  load  resistance  (Qm  for  the  ceramic  decreases  with 
increased  drive  level). 
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Figure  A-3.  Quality  Factor  versus  Amplifier  Impedance 
(f|  = 20  kHz,  Qfo  = 500,  k" = 0.3,  and  Cb  = 10'9  F) 

The  results  of  the  analysis  can  only  be  used  in  a qualitative 
sense  because  the  constraints  imposed  on  the  amplifier  impedance  did 
not  include  any  frequency  dependence  nor  the  fact  that  the  amplifier 
was  actually  used  in  a pulse  mode.  A more  rigorous  analysis  could  have 
accounted  for  the  reactive  components  of  the  amplifier  impedance,  but 
quantitative  results  could  not  have  been  established  without  an  exten- 
sive measurement  program  to  determine  the  true  amplifier  impedance  com- 
ponents at  different  frequencies  and  load  levels.  Since  this  was  not 
done  in  the  present  circumstance,  it  will  be  assumed  that  the  influence 
of  the  amplifier  impedance  on  determining  the  quality  factor  of  the  ce- 
ramic rings  at  resonance  is  negligible. 
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Table  A-l.  Measured  Values  of  Auplifier  Impedance  (200-W  Macintosh) 


Frequency 

f 

Open-Circuit 

Voltage 

Load 

Resistance 

Voltage 

Across 

Plate 

Current 

Anplifier 

Inpedance 

r 

(kHz) 

%c(t) 

(Vrms) 

rl 

m 

|e(t)| 

(Vrms) 

(S) 

(!!i  , 

o. 5 ki> 

i 

2 

450 

91.5 

7 n 

185 

51 

90 

450 

51.5 

i.  u 

100 

140 

400 

32 

500 

180 

110 

14 

1 

2.05 

450 

91.7 

10. 

190 

51 

88.5 

450 

58.5 

100 

137. 

375 

40 

500 

185. 

160 

13.5 

1 

2.15 

450 

75.7 

18.5 

165 

51 

100 

88 

137 

450 

400 

44.6 

20 

SOO 

165 

210 

— 

THE  CURRENT -DECAY  CASE 

The  influence  of  the  amplifier  impedance,  under  the  same  restric- 
tions introduced  in  the  previous  case  (namely,  that  the  impedance  could 
be  represented  by  a real  lumber- Ra),  was  also  analyzed  by  using  Laplace 
transformation  techniques.  The  circuit  in  figure  A-l  still  applies 
with,  as  before,  Za  replaced  by  Ra  and  eoc(t)  set  equal  to  zero.  The 
two  loop  equations  for  the  circuit  then  become  simply 

0 " ilRa  + ~T  [/ildt  ' /i2dt]  (A -11) 

b 

and 

dd2  1 1 

0 - + Ly  /i2dt  + -±-  [/i2dt  - / ijdt],  (A -12) 

i b 

where  ij  and  are  the  instantaneous  currents  in  the  two  loops.  The 
intent  is  to  determine  the  decay  behavior  of  ij,  which  was  monitored  as 
the  voltage  drop  across  a 1-0  resistor  placed  in  series  with  the  ring 
impedance  and  the  low  (ground)  terminal  of  the  amplifier.  This  addi- 
tional resistor  appears  in  loop  one  and  may  be  effectively  accounted  for 
by  letting  Ra  * (Ra  ♦ l)  in  any  calculations. 
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The  result  of  applying  the  Laplace  transform  to  equations  (A-ll) 
and  (A- 12)  and  rearranging  terms  is 

‘l(S)  {A) 

■ ^ [ij'1’  (°’>  - *{'I)  (»*)]  CA-i: 

and 

i2(s)k*  sly* 

* hV0*’  - ^ (k1^0*)) 

- [1!'1)',‘)  - k'V’]  • »-i‘ 

Initial  conditions  must  be  specified  for  the  circuit  and  will  be 
established  by  using  the  information  contained  in  figure  A-4.  Since 
e(t)  = E0cos(u)^t)  appears  across  both  branches  of  the  circuit  during 
steady-state  operation,  then 


i 2 1 c ) = -jp-  cos(^t)  at  = 


The  voltage  across  the  inductance  Ly  is 


E,  o . , E., 
wrLY  “R^- 


which  is  equal  and  opposite  in  sign  to  the  voltage  across  the  capaci- 
tance element  Cy.  Therefore,  at  t = 0,  the  intial  conditions  can  be 
expressed  as 


i^n(0+)  = qY(04)  = 0, 

V0+>  - ~rr  • 


(A-17) 


F" 
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and 


i{_1)(0+)  - i^1}(0+)  = qj(0+)  - q2(C+)  = qb  = , 

where  q designates  electrical  charge. 


t * o 


Equations  (A - 1 3 ) and  (A-14)  can  now  be  solved  for  Ii(s),  by  using 
the  initial  conditions  in  equation  (A-17),  to  obtain 


I^s)  = 


where 


o 

tRY 


s ♦ 


J^r_ 


Is2  * A. 
W \ % 


s + 0) 


, (A-18) 


/ 1 \(  2 “r  I2\ 

{s'~\  s'"r 


r vrn? 
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and 


1 


y 


i - R C,  . 
a b 


The  current  ij(t)  is  found  by  taking  the  inverse  Laplace  trans- 
form of  equation  (A-18),  which  has  the  form 
E 


ix(t) 


-^|-|2eat  [AlCos(bt)  - A2sin(bt)]  - A3eCt|,  (A-19) 


where  a,  b,  and  c are  numbers  associated  with  the  poles  of  equation 
(A-18)  and  where  Aj,  A2>  and  A3  are  coefficients  determined  by  using 
residue  theory  at  the  poles.  The  pole  values  of  equation  (A-18)  were 
determined  for  different  values  of  the  time  constant  t by  using  the 
polynomial  root  finder  program  on  a desk  calculator*  for  the  same  val- 
ues of  fj,  Q^,  k,  and  Cj,  as  were  used  in  the  first  case.  The 
resulting  values  were  fed  into  another  program  to  determine  the  coeffi- 
cients Aj,  A2,  and  A3.  All  values,  in  turn,  were  incorporated  into  a 
plotting  routine  program  for  equation  (A-19),  and  the  effective  quality 
factor  Qm  was  determined  from  the  individual  decay  curves  for  Ry  = 100  SI 
and  E0  = 10  V.  In  all  cases  the  last  term  in  equation  (A-19)  was  found 
to  be  of  minor  significance,  and  Qm  was  actually  determined  from  the 
leading  term  in  the  equation.  Table  A-2  lists  some  of  the  values  ob- 
tained by  using  this  approach,  and  it  should  be  noted  that  the  frequency 
coefficient  lies  between  wf  (1.2566  x 105  Hz)  and  io*  (1.3173  x 10=  Hz) 
as  Ra  increases.  The  Qm  values  obtained  were  superimposed  on  figure  A-3 
to  demonstrate  the  agreement  between  the  two  methods. 

The  agreement  between  methods  should  be  as  good  in  general  as  that 
shown  here  since  the  same  circuit  and  circuit  elements  were  analyzed  in 
both  cases.  Even  though  the  behavior  of  Qm  versus  amplifier  impedance 
is  highly  nonlinear  over  the  broad  range  of  values  considered,  at  any 
discrete  value  the  circuit  elements  are  fixed  and  linear  analysis  is 
valid.  The  agreement  demonstrated  that  current  decay  measurements  could 
indeed  be  used  to  predict  the  quality  factor  variation  for  ceramic  rings 
as  long  as  other  practical  considerations  mentioned  in  the  text  (such 
as  dependence  of  Qm  on  decaying  stress  level)  are  not  disregarded. 


•Hewlett-Packard  model  9820-A 
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ft 
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Table  A-2.  Values  for  versus  R 


Amplifier 

Impedance 

R 

a 

Decay 

Coefficient 

a(xl02) 

Frequency 

Coefficient 

b(xl05) 

c 

Effective 
Quality  Factor 

Qm 

10"4 

-1.25 

1.2566 

-1013 

500 

10'1 

-1.26 

1.2568 

O 

H 

o 

H 

1 

499.8 

25 

-1.45 

1.2569 

-4xl07 

433 

50 

-1.65 

1.2569 

r^ 

o 

H 

X 

rs* 

380 

102 

-2.00 

1.2569 

-lxlO7 

308 

S.OxlO2 

-5.00 

1.2570 

-2xl06 

126 

103 

-9.00 

1.2582 

-lxlO6 

70 

8. OxlO3 

-32.31 

1.2890 

-12xl04 

19.4 

105 

-5.70 

1.3170 

-9xl03 

110 

io7 

-1.30 

1.3171 

-9X101 

507 

1010 

-1.25 

1.3173 

-9xl0-2 

524 

= 2n(20xl03)  Hz 

P-Y  = 100 

n 

Qr  = 500 

E = 10  V 
0 

k * 0.3 

Cb  = 10 -9 

F 

A- 11/A- 12 
Reverse  Blank 


I 
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APPENDIX  B 


DETERMINATION  OF  TOTAL  CURRENT 


The  expression  for  the  total  current  I into  the  ceramic  ring  is 
given  by  equation  (9)  in  the  text,  which  is  (repeated  for  convenience) 


I = jcoA 


d31d31 


r-  . ; »_T  / , d31d31 

E3  Ja)A  33  1 1 T E I E3  ’ 
\ e33Sll  / 


(B-l) 


If  it  is  assumed  that  1/Qm  < 1 or  that  sfi  can  be  represented  as  a real 
number,  then  at  resonance  equation  (B-l)  reduces  to 


“rAd31d3lQM 

E 


T i jx  n n . E.  s _ 

1 = E E*  + JU)rAe33E3* 


(B-2) 


where  the  blocked  permittivity  is  defined  for  linear  conditions  as 


s T ri  ,2, 

33  ' e33  " k J ' 


(B-3) 


The  quantities  in  equation  (B-2)  can  be  further  redefined  to  yield 


1 ■ -sr  * j“rCbE0’ 


(B  —4 ) 


which  shows  that  the  input  current  is  composed  of  a blocked  component 
(lb  = jwfCbEo)  and  a motional  component  (Imot  = eo/rY)»  »*ere  E0  is  the 
magnitude  of  the  driving  voltage  during  steady-state  operation.  This 
combination  is  easily  represented  by  the  two-branch  circuit  in  figure 
A-l  (in  appendix  A)  to  the  right  of  terminals  a and  b.  Since  the  com- 
ponents of  the  input  current  are  in  quadrature  at  resonance,  and  since 
the  magnitude  of  the  input  current  was  the  quantity  measured  during  the 
experimental  phase  of  the  program,  the  relationships  depicted  in  figure 
B-l  will  be  used  to  determine  the  motional  component  of  the  current. 
From  figure  B-l  it  is  obvious  that 


B-l 
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and  that 


1/2 


(B-5) 


tan  0 


(B-6) 


F 

Figure  B-l.  Components  of  the  Input  Current  at  f = f 

Now,  the  blocked  component  of  the  input  current  varies  directly  with  the 
applied  voltage,  equation  (B-4)  (that  is,  |l^|  = wpCFU  - k2)ED>*  where 
CF  is  the  low  frequency  capacitance  of  the  rings).  Therefore,  the  max- 
imum correction  would  be  associated  with  the  highest  field  level  used. 
For  ring  A at  a field  of  E3  = 12.5  kV/m, 

Ife  = 211(18. 374  x 103)  (890  x 10‘12) ( 125)  (1 . 27)  = 16.3  mA.  (B-7) 

The  measured  input  current  at  the  same  field  level  is  I = 239  mA,  which 
yields  a ratio  of  |lb|/|l|  = 0.068,  or  | Imot j = (0.998) | I | ; this  value 
indicates  that  the  blocked  component  of  the  input  current  may  be  ne- 
glected for  all  practical  purposes  in  this  report. 


"Although  CF  is  known  to  increase  with  increased  electric  field  level, 
the  field  levels  during  these  experiments  were  considered  to  be  suffi- 
ciently low  that  the  capacitance  remained  essentially  constant. 
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